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Interpretation Geometrique du Probleme de 
Compensation Non-Lineaire 

Abdelmajid Ben Hadj Salem 



1. Introduction 

Dans un article |1J E. Grafarend et B. Schaffrin ont etudie la geometrie 
de la compensation ou I'ajustement non-lineaire et ont presente le cas du 
probleme d'intersection plane en utilisant le modele de Gauss Markov, par les 
moindres carres. Le present papier presente les principes de la geometrie de la 
compensation non-lineaire par la methode des moindres carres en s'appuyant 
sur le Lemme de Pazman [2]. 



2. La Geometrie Non Lineaire du Modele de Gauss-Markov 

Le modele non lineaire de Gauss-Markov est defini par : 
(1) C{X) = L-e; eeM{0,T) 

avec : 

- L : le vecteur des observations (n x 1) = (Li, L2, Ln)"^ , 

- X : le vecteur des inconnues (m x 1) = [Xi,X2, ..^Xm)^ , 

- e : le vecteur des erreurs (n x 1) = (ei, 62, .., e„)"^ suit la loi normale 
7\A(0, r) avec E{e) = et L = E{ee^) la matrice de dispersion ou variance, 
on prendra T = cTq.P"^. P est la matrice des poids et (Jq une constante 
positive. 

- C, : est une fonction donnee injective d'un ouvert U C — ?• M" et m < n. 



Remarque : dans le cas d'un modele lineaire, la fonction C, = A.X ou A 
est une matrice n x m. 

On note Im( = {(^{X)/X G U} I'image de U par la fonction ImC, est 
une variete de dimension m verifiant les conditions : 

(1) : les vecteurs — — ,— — , 7; sont Imeairement mdependants en 

dXi 0X2 oXm 
chaque point X , 

(ii) : les fonctions — — sont continues sur U pour i, j € |1,2, ...,m|. 

^ ' dXidXj ^ ,J \ , , ^ f 



3 



On introduit un produit scalaire : 

(2) <Ci,c2>=cT-p-C2 

D'ou la norme du vecteur = {(^i,C,2, Cn)^ • 



(3) \\cf=<cx>=C^.P■C = ^P^.C 



2 

i 

i=l 



dans I'espace vectoriel M" en prenant la matrice de poids P une matrice 
diagonale . 

Alors la solution par les moindres carres X sera definie par : 

(4) \\L - aX)\\ = min{\\L - C{X)\\ / X G U} 

Cette condition est exprimee par les equations de Lagrange-Euler soit : 

(5) A_\\L - C{X)f = pouriG{l,2,...,m} 
En effet, on veut minimiser la fonction : 

(6) F{X) = F{Xi,X2,..,Xm) = \\L - C{X)\\ = \\L - ({Xi, X2, X^)\\ 

Comme F est une fonction positive, minimiser F c'est aussi minimiser F^, 
soit J{X) = F'^(X). En appliquant les equations de Lagrange-Euler, on 
obtient : 

dJjX) _ ^ dJjX) _ 
dXi dXi 

soit : 

(7) ^\\L-C{Xi,X2,...,Xm)f = pouriG{l,2,...,m} 
or : 

\\L - C{Xi,X2, = {L- C{Xl,X2, ...,Xm) f.P.{L - C{Xi,X2,...,X^)) = 

(8) C{Xf.P.C{X) - 2L^.P.C{X) + L^.P.L 
Soit : 

ou encore : 

^ = 2(C(X) -L)^.P.^ pour. E {1,2,...,^ 

ce qui donne en utilisant ([5]) : 

(10) <L-C(X),^|^ >=0 pouriG{l,2,...,m} 
. ^ dCiX) 

(11) ou < e, >= Q pour i G {1,2, ...,m} 

oXi 
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Geometriquement, cela veut dire que le vecteur erreur e = L — C{X) est 
perpendiculaire (produit scalaire nul) au plan tangent de la variete ImC, au 
point C,{X) (s'il existe). 

Pour le cas non-lineaire, la condition (|10p est necessaire mais non suf- 

fisante. Pour obtenir le minimum, il faut que la matrice (— — ),«,?' G 

^dXidXj 

{1,2, ...,m} soit definie positive. 



3. Interpretation Geometrique 

Dans cette section, on va clarifier I'interpretation geometrique de fagon 
que la solution de (|1U|) soit localement unique. 



Considerons la matrice m x m definie par : 



(^9)r(x)-(r ) ..err _^ 9C{X) ddx) (a = i,2,..,m\ 

(12) G{X)-{G^,) avecG.,-<^^,^> ^ 2, m j 
Or : 

(13) ds^ = Go,pdXc,dXp 

represente la metrique de la variete ImQ. La matrice G{X) = {Gap) est 
appelee en terme statistique la matrice d'information de Fischer. 



Introduisons la matrice B definie par : 
(14) 

B(X,L) = (B„,) avecB„» = l^||L-((X)f (^rj;,^;:;:) 

or : 

\\L - aX)\? = {L- aX))^P{L - C(X)) = (L^ - e{X)){PL - PC{X)) = 
(15) L^PL - 2L^PC{X) + (^{X)PC{X) 

D'ou : 

^ -{L^PL - 2L^PaX) + C^(X)PC(^)) = -21^ P^ + 2C^{X)P^^^^^ 



a 



dXa ^ V y QJ^^ ^ V / QJ^^ 

(16) -2(L^-C^(X))P^ 
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Done : 



2dXadX^" dX^dXp dXp dX^ ■ ^ V / dX^dXp 



dXp dX^ ' ^ V dx^dXp 

im < ^^^^^ > < L ax-) ^'^(^^ > 

Soit : 

(18) B^, = G^,-<L-aX),§^^> 
Posons : 

(19) H(L,x)=(K,)=«L-axi§^^» 

c'cst-a-dire : 

(20) B = G-H 

Revenons a Im({X) = {C(^) /X G U}. Soit une ligne geodesique de ImC,{X) 
passant par un point C, = C,{X) parametree par son abcisse curviligne s, on 
a alors : 

(21) x{s) = aX{s)), sG[si,S2] 

ou X{s) decrit une certaine courbe dans le domaine U C W^. 

Le vecteur : 

(22) x'W = "-^ 

represente le vecteur tangent a la ligne geodesique au point C,{X{s)) de 
Im({X). Ce vecteur verifie : 

(23) llx'(^)f = l 

Par suite, la derivee de ce vecteur par rapport a s est un vecteur orthogonal 
a x'(s) done orthogonal a Im({X{s)) au point ({X{s)) : 

(24) ^^s) = ^±x'is) 

c'est-a-dire parallele au vecteur normal a la surface ou la variete ImC,{X) et 
on retrouve la propriete que x(s) est une geodesique. 

Remarquons que pour une ligne geodesique, la courbure geodesique est 
nulle et la courbure normale coincide avec la courbure de la courbe %(*) 
soit : 

(25) p{s) ^ 



llx"(^)ll 
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le rayon de courbure. Appelons : 

(26) n{s) = ^^^=r{s).p{s) 

c'est un vecteur unitaire perpendiculaire au plan tangent a la surface ImC,{X). 

D'apres I'equation (llOp . au point C(-^)) 1^ vecteur e est perpendiculaire a 
L — C{^)- Appelons alors : 

(27) i^(CW) = \Z/ Z e M"avec < ^, >= a = l,2,...,m 

On done e = L — C(^) £ K. Ce dernier est un espace vectoriel de dimension 
n — m orthogonal a ImC,{X) au point C(^)- On ^ aussi x"(s) G K. 

3.1. Lemme de Pazman. On pent maintenant enoncer le lemme de Paz- 
man (1984, [2j) comme suit : 

Lemme de Pazman : Pour tout vecteur d' observation L G , et toute 
solution appropriee X des equations : 

<L-C{X),^^>=0, a = l,2,...,m 



les conditions suivantes sent equivalentes : 

1)- La matnce B{X,L) = G{X) - {< L - C(^),5f^ 
positive. 

2)- Pour toute ligne geodesique xiX{s)) verifiant : 

= Q{X{s)) 

On a I'inegalite : 

(28) <L-Q{X),n{s) > < p{s) 



A. La condition 1 : 

En effet, supposons que la matrice B{X,L) est definie positive c'est-a- 
dire : 

(29) VF G , y / ^ Y'^.B.Y > 
Prenons alors .Y = x'i^)- On a : 

(30) x'^{s).B{L,X)-x'{s)> 
Comme B = G{X) - H{L,X), on obtient : 

x'^mGiX)-HiL,X)).x\s)> 

soit : 

(31) x'^{s).G{X)x'i-s) - x'^miL,X).xis) > 



Or pour s : 



i=i * i=i * 

en notant ~ ^^'j'^^ les composantes de 'x'{s) dans le plan tangent a 

ImQ au point C,[X{s)). Comme : 

i=i * j=i * 

(33) 2^ ( < — — , > I Xj(s) = X (s) -G-x (s) = 1 

i,j=l * J / 

Prenons s = s, alors (j3ip deviant : 

(34) x'^{s)H{L,X).x'{s)< 1 

Comme B = G — H done la matrice H est exprimee dans la base de B soit 
( En utilisant ([32]), le nombre reel x'^(s)i?(L, X).x'(s) s'ecrit : 



V 

(35) 

i=m I j=m 

x'^miL,X).x'{s) = ^x-(^-) E ^^^-^K^") 1 = E X[{s).x',{s).hj 

i=l \ j=l J i,j=l 



On remplace hij par 



Par un calcul simple, I'equation (I35p devient : 

(36) x'^(^-)i^(i,^).x'(^-) =< i - C(^), E > 

Maintenant, on va s'interesser au membre droit du produit scalaire de I'equa- 
tion (I36p . En differentiant I'equation (I32p par rapport a s, on obtient : 

i i j 

Alors on a pour s = s : 

<.-c,.),."(.)>^<.-ax,,E^™.E.:«E»f^> 

i i j 

(38) 

E^.<-c(.,™>.Ex:(.E<-c,n».^ 

7 1 l' 



Or en utilisant I'equation (llOp . le premier terme de la deuxieme ligne de 
I'equation precedente est nul : 

i 

et comme : 

X'(«) = ^ 
ds 

Alors I'equation (j38p devient : 

<L- ann-s) >= E x^{s) E < ^ - c(^), ^m^-^K-) >= 

(39) ^^,;(,),;,.)<L-cm.^^> 

Or le deuxieme membre n'est autre que I'equation (I36p . En utilisant (I34p . on 
obtient : 

(40) <L-CiX),xHs)><l 
Or : 

n(s) = />(s).x"(s) 

D'ou : 

(41) <L-C{X),n{s)> < p{s) 
B. La condition 2 : 

On suppose que reciproquement, on a pour toute geodesique x(s) = 
C(^(s)) de ImC passant par le point x(^) = C{X{s)) verifiant : 

(42) <L-C{X),n{s)> < p{s) 
et : 

<L-C{X),^^>=0 pouriG{l,2,...,m} 



lx'(^-)f = 1 



^(^) = ^^=P(^)X"(«) 



telque x'i^) verifiant : 
Comme : 



le remplagant dans I'equation (142 p . on obtient : 

(43) <L-aX),p{s)xHs)> < p{s) 
et en simplifiant par p ^ 0, soit : 

(44) <L-C(X),x"(s) > < 1 



Comme : 



1=1, m 

dc 



i=l.m 



D'ou en derivant une deuxieme fois par rapport a s : 

- V MM ^C(^(^)) , V v' r \- d\{x{s)) dXj_ 
^^^"^ ds ■ dx, ^^'^^ ^ dx,dXj ■ ds 

i i 3 

En remplagant x"('5) dans (j44p . on obtient : 
(46) 

i * J 

ou encore : 
(47) 

i i j I- J 

Or la premiere somme est nulle en vertu de I'equation (llOp . II reste : 

(48) Y.Y.^^{s)m < L - axi %%^^ > <1 



Comme : 



f( N ^ ^x(g) ^ (9C rf^i(s) 

«=l,m 



1=1, m 

et x'(s) verifie ||x'('S)ll = 1 car c'est un vecteur miitaire tangent a la geode- 
sique xi^)- Done 



(49) < x'{s),x'{s) >=1^< Y E 



— dXi' ^ 'dx. 

i=l,m j=l,m ■' 



>= 1 



soit : 
(50) 



E E xK-) < ^ > x'M = E E xK^)G.,x;.(s) = 1 

i=l,mj=l,m * i=l,mj=l,in 
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et ce-ci n'est autre que : 

(51) x'^{s).G.x'is) = 1 

En utilisant I'equation (|48p . on a : 

i j « J 

ou encore : 

(53) < x'^is).G.x'{s) -Y^Y.^^(s).h,j.X-is) 

i j 

Finalement, nous obtenons : 

(54) < x"^{s).G.x'{s) - x'^{s).H.x'{s) 
C'est-a-dire pour tout vecteur Y = x'i^) 7^ du plan tangent de Jm^ : 

(55) Y^ .B.Y > =^ la matrice B est definie positive 

C.Q.F.D 



Maintenant, on peut dire quand X solution de (llOp est solution des moindres 
Carres en reference a I'equation ([4]) et ce a partir du corollaire suivant : 

Corollaire : Si X est solution de : 

(56) <L-C{X),^^>=0 pouriG{l,2,...,m} 

avec : 

(57) ||L-C(^||<r 

oil r designe le rayon de courhure minimum de la variete Im( defini par : 
(58) 

r = inf {py.{s)/xis) = C(^(s)) une geodesique passant par ({X),\/X € U} 
Alors X coincide avec la solution des moindres carres X = X{L). 

Demonstration : 

Comme les deux vecteurs e = L — C(^) et n{s) sont ortliogonaux a ImC, 
lis sont colineaires et comme n{s) est un vecteur unitaire alors leur produit 
scalaire < L — (^(X),n(s) > est plus petit ou egal a \\L — C(^||- Or ce 
terme est plus petit que r d'apres (I57p . et comme r est le plus petit rayon 
de courbure, on a r < p^is) pour toute geodesique x(s) passant par C(^)- 
Ce-ci est traduit par I'equation : 

(59) <L- C{X), n{s) >< \\L -C{X\\<r< p^{s) 
De cette derniere equation, on a : 

(60) <L-aX),n{s)><Px{^) 
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En utilisant le lemme de Pazman cite ci-dessus, la matrice B est definie posi- 
tive done C(^) 6st un minimum strict [3]. Or on a suppose que I'application 
est injective (si Ci^i) = ==^ = ^2), alors X coincide avec la 

solution des moindres carres X = X(L). 
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